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ABSTRACT 

In this paper, we revisit in detail the effects of primordial magnetic fields on big 
bang nucleosynthesis (BBN) including a discussion of the magnetic field geometry and 
the anomalous magnetic moment. The presence of magnetic fields affects BBN by (1) 
increasing the weak reaction rates; (2) increasing the electron density due to changes to 
the electron phase space; and (3) by increasing the expansion rate of the universe, due 
both to the magnetic field energy density and to the modified electron energy density. Of 
the effects considered, the increase in the expansion rate due to the magnetic field energy 
is the most significant for the interests of BBN. The allowed magnetic field intensity at 
the end of nucleosynthesis (0.01 MeV) is about 2 x 10 9 G and corresponds to an upper 
limit on the magnetic field energy density of about 28% of the neutrino energy density 
(p B < 0.28p„). 

PACS number (s): 98.80.Cq, 98.62.En 

1. Introduction 

Big bang nucleosynthesis (BBN) provides an unique quantitative window for processes 
occurring in the early universe^ between temperatures of lMeV and O.OlMeV. The agree- 
ment between the light element abundances predicted by BBN and observations strongly 
constrains the dynamics of the universe at this epoch including the presence of strong mag- 
netic fields. A primeval magnetic field existing during nucleosynthesis would have three 
major effects on BBN: (i) it would alter the weak interaction rates, (ii) it would modify 
the electron densities in phase space, and (iii) it would increase the cosmological expansion 
rate. Some of these effects were examined by a number of authors^ and most recently 
by Cheng, Schramm, and Trurant 3 ! , Grasso and Rubinstein^ , and Kernan et al^ . In this 
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paper, we revisit our earlier analysis ^ and find reasonable agreement between subsequent 
work by different authors! 4 ' 5 ] and our present results. Although some slight differences 
remain, the basic conclusions seem unambiguous. Here we also show that the effects of 
the spatial distribution of the magnetic fields and the anomalous magnetic moment do not 
affect significantly the results. 

2. Three major effects of B fields on BBN 

Jedamzik et ah 6 ! have shown that neutrino decoupling effectively damps all magne- 
tohydrodynamic modes up to the scales around a tenth of the Hubble radius at neutrino 
decoupling. If Id is the largest scale over which the magnetic field becomes spatially ho- 
mogeneous due to neutrino damping, then Id ~ O.lif -1 at T ~ 1 MeV. This implies that 
if there are magnetic fields present during BBN their spatial distribution is very smooth 
on scales smaller than Id and the field can be taken as constant within these scales. 

The magnetic field spatial distribution needs to be taken into account if Id is smaller 
than the length scales over which reactions and mixing occur during BBN. The relevant 
scale to be compared to Id here is the largest mixing length which corresponds to the 
neutron diffusion length, d n . Jedamzik and Fuller' 7 ] showed that d n (l MeV)<l m while 
the horizon MeV) ~ 10 8 m. Since Id ^> d n , the magnetic field is constant within 

correlated volumes and will be taken as constant below. 1 We also assume that the field 
is randomly oriented within each volume of radius Ib such that the expansion rate is not 
anisotropic and that Robertson- Walker metric is valid. 

In an uniform magnetic field with magnitude B chosen to lie along a z-axis, the 
dispersion relation for an electron propagating through the field is 

E = \p\ + m\ + 2eBn s ]^ + m e K, (2.1) 

where n s = n + \ — s z , (n s = 0, 1, ...), n is the principal quantum number of the Landau 
level and s z = ±1/2 are spins, e is the electron charge, p z the electron momentum, m e the 
rest mass of the electron, and k is the anomalous magnetic moment term! 9 ! for an electron 
in the ground state (n = 0, s z = 1/2). For relatively weak fields (i. e., S<7.575 x 10 16 

1 Recently, Grasso and Rubinstein! 8 ! assumed that the magnetic field during BBN had 
fluctuations on the scale of the horizon at the electroweak transition which is of the order 
of the diffusion length at the end of BBN. If damping due to neutrino decoupling were 
not as effective as found in Ref. [7] , the spatial variations of the magnetic field could affect 
BBN outcome. 
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G), k = — f^^fl, while for stronger fields, k = f^ln 2 ^? ) 2 , where a e = The number 
density of states in the interval dp z for any given value of n s in the presence of magnetic 
field is described bj^ 

(2-*».o)^dp,. (2-2) 

We now discuss how Eqs.(2.1) and (2.2) affect BBN in detail. 

2.1 Weak reaction rates: 

The weak interaction rates in a constant magnetic field without QED correction have 
been derived by Cheng, Schramm, and Truranl 11 !. In an expanding universe, the B field 
evolves as -R -2 , where R is the scale factor of the universe. During BBN, R oc T" 1 , 
where T v is the neutrino temperature. Let Bi represent the magnetic field at an initial 
temperature = 1 MeV, 7; = Bi/B c and 7 = B/B c , where B c = m 2 e /e = 4.4 x 10 13 G is 
the critical field at which quantized cyclotron states begin to exist. The magnetic field at 
any temperature T can then be written as 

B = B cl A 2 or 7 = 7,(^) 2 - (2-3) 

With this notation, the rate for the reaction n + e + — > p + P e , is given by 
+ 3a 2 )mh*T 2 ^ ro . .f°° de(e - «) 



_ ^(l + 3^)m^ ^ / 

* n s =0 



' v / T+2^ir+ K a/(c - k) 2 - (1 + 27n s ) 



1 ( e + g )2 e ( e +flf)z„+^ 

(1 _|_ e eZ e +4> e ^ (I + e (e+q)Z v +4, v -j ' 



(2.4) 



where ^(1 + 3a 2 )m^/2yr 3 ~ 6.515 x 10 _4 sec _1 , g v = 1.4146 x 10" 49 erg cm 3 , and a 
g A /g v ~ -1.262J 12 ] 

Similarly, for the reaction n + v — > p + e~, the rate is 



^(l + 3a 2 )m^T 2 ^ f9 , , r de( 



■ ? n s =0 



V i +27ns+K y/(e- n) 2 - (l + 27n s ) 



1 (e - g)V Ze+ * e 

X (1 _|_ e eZ e +<£ e ) ^ ZjZ e (e-g)Z„-<^ 

^(l + 3a 2 )m^T 2n v^ x ro , , T cfc(e - «) 



' l 'S max /»(^ 

A/ 



(2.5) 



^T 2 ^ 1 n ^ 7^1+2^7+. v / (^-^) 2 -(l + 27^) 

1 (e - g)V Ze+ * e 



x 



(1 _|_ e eZ e +<£ e ) ^ _|_ e (e-q)Z v -4>„y 
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and for the reaction n — > p + e + z/ e , we have 

_#(l + 3aVM V*"ro * i f cfc(c - «) 



' < s max pq 



1 (g - e ) 2 e eZ e + ^»e 

(1 _|_ e eZ e +(j} e '-j ^ _|_ e (q—e)Z l/ +4> l/ ^ ' 

The total weak reaction rates for the conversion of neutrons to protons is simply the sum 
of the above rates 

! n s =0 ^Vl+^n.+it ^ 

1 l t + q f e {e+q)Z v+ d, v ( e _ q j2 e eZ e +<t> 

X [( e -«)2 - (l + 2 7 7i s )]5 ^ l + e (9+ £ )^+^ + l + e (e-g)z„-*J ' 

The parameters used above are defined as e = — , q = m "~ mp , Z e = Z v = 
4> e = iff, and <f> u = %j¥-, where m n and m p are the rest masses of the neutron and 
proton respectively, T e is the temperature of the electrons, = e, v) is the chemical 
potential of the electron or neutrino, and <f>i(i = e, v) is the degeneracy parameter. 

For /3-decay processes to occur in the presence of a magnetic field, the quantum number 
n s has to satisfy 



a/1 + 2^n s + k <q, or n s < n smax = Int 



(2.8) 
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where n smax is the largest integer in — ^ ~ 1 . (Note that our expressions for the total 
rate (2.7) here is different from that given in Ref [4].) 

The inverse total reaction rate of the conversion of protons to neutrons is computed 

as 

A =p-i z e\ 

In order to better elucidate the effects of the field on the rates, we calculate analytically 
the variations of the reaction rates with respect to changes in the field, and we obtain 

lim -± oc E t 2 " *n s o] / dfc/+(e)« + ( 7 , e) > 0, (2.9) 

7 \ pOO smax pq 

lh ?o^T* Et 2 -^] / */-(c)«-(7,c)- E [2- 5 «.o] / dfe/_(e)«_( 7 ,e)>0, 

(2.10) 
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n s =0 

and 



7 \ ^smax pq 

lim-^oc ^[2-5 ns0 ] / dkf_(e)u_( 1 ,e)>0, (2.11) 



J\ 00 /-oo 

lim-p^oc ^[2-5 ns0 ] / dfc[/ + (e)« + (7,e) + /_(e)«_( 7 ,c)]>0, (2.12) 



n s =0 

where 

2 

e= (k 2 + l + 2 1 n s ) 1 / 2 + K, fc = % 



(1 _|_ giCeZe+^e)^! _|_ e -[{q±t)Z v +4> v \y 



W±l7,Cj e(e±(/) T e 1 + e ±(^«+*e) e(l + e -[(«±^+<M) ' 

We also computed Eqs.(2.9 -2.12) numerically, for various of 7$ and T. Both calcu- 
lations show that, independent of the temperature, the presence of a magnetic field does 
increase all the weak reaction rates, including the total neutron depletion rate. This result 
is consistent with the findings in our previous works and with Grasso and Rubinstein's 
recent calculations, M but inconsistent with Kernan's recent statement^ that the rates of 
2-2 processes decrease as the field increases. At very high temperatures T 3> 2.5 MeV, 
such effects are insignificant because the inverse reaction rates also increase with the field 
and are not much suppressed by the factor exp(— qZ e ). When the temperature drops to 
a point where the reactions n + e + — > p + V and n + v — > p + e~ begin to freeze out and 
the neutron /3-decay process dominates, then the total rate increases with the magnetic 
field. However, if the primeval field is not strong enough to begin with, then as the uni- 
verse expands, it becomes too weak to affect the reaction rates at low temperatures. Our 
numerical calculations reveal that for the magnetic field to have significant impact on the 
reaction rates, 7^10 3 or Bi>AA x 10 16 G at 1 MeV. As we discuss below, the effect due 
to the change in expansion rate is already significant for 7^ 10 and dominates over the 
change on the reaction rates. 

2.2 Electron density phase space 

In a magnetic field, the phase space and energy density of electrons are modified. The 
number density and energy density of electrons (n e and p e ) over phase space as a function 
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of magnetic field strength are given by 

,3„, t2 00 



n e = 2 me T5; Y 0) r cfe - (e ^ (2-13) 



and 



(2.14) 

Correspondingly, the pressure of electrons is 



p _ 2 ™^ ^ /~ ^ - «)* - (1 + 2-yn.) 

p « - 2 W ,^ (2 " *"■») ' (2 ' 16) 

These expressions will reduce 



3 r r eV^i 



n e = I de- 
11 2 J 1 1 



_|_ e eZ e +(f) e ' 

1 <- x -v^ r 



(2.16a) 



p e = ^ I de i v \ , : , (2.166) 



and 

P = m e r ^ (^ 2 - !) 3/2 



/ de— (2.16c) 

if the magnetic field is absent. 

The dependences of n e , p e , and P e on the magnetic field can be seen analytically to 

be 

lim^ocf (2- 5 n 0) r dk V^H 1 " ~, ^fy ^ d > °> ( 2 -17) 

lim^oc V(2- 5 ns0 ) 1^ dk- "-^-—Al + ^r- „ in f 7 \ s ^ ] > 0, (2.18) 

7^0 d'j ^ 3 ; J (l + e eZe +<^) L e 2 e(l + e"( eZe +^))) 

and 

lim^oc y (2-8 n 0) f°d£;- ^ — - [1 - ^ - — ] > 0. (2.19) 

7-0 d 7 ^ v nsU Vo e(l + e eZ =+^) e 2 e(l + e -( eZ -+^)) J V ; 

where we assumed non-degenerate neutrinos (0„ = 0). These expressions indicate that, 
in the presence of magnetic fields, due to the large Landau excitation energy and the 
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decreased cross-sectional area of each Landau level, all of the electron thermodynamic 
quantities, such as n e , p e , andP e , increase with increasing field strengths. This, in turn, 
causes a decrease in all of the weak interaction rates and changes the temperature-time 
relationship in BBN calculations. Furthermore, it results in an increase in the final neutron 
to proton ratio and the abundances of the light elements. Such an effect becomes significant 
for 7i>10 3 and is sub-dominant to the effect discussed below. 

2.3 Effects on the Expansion Rate 

The expansion rate of our universe is given by 



1 dR 8itG 

^iW— » < 2 - 20 > 

where G is the gravitational constant and p is the total energy density. This can be 
expressed as p = p 1 + p e + p v + p b + p B , where p e = p e - + p e +, p v = p Ue + p Vfi + p Vr + 
Pv a + + Pv T > an d the subscripts 7, e, v e , z/ M , z/ r , 6, and B stand for photons, electrons, 
e-neutrinos, /U-neutrinos, r-neutrinos, baryons, and magnetic fields. 

The presence of magnetic fields alters the expansion rate by the added energy density 
of the magnetic field 



and the change in the electron energy density which we can write as 

p e = p e (B = 0) + 5p e (2.22) 

During nucleosynthesis, e + e~ annihilation transfers entropy to the photons but not 
to the decoupled neutrinos. The neutrino temperature then follows T v oc -R -1 , i.e., 

dT, 
dt 

while the photon temperature satisfies 



= -HT V (2.23) 



dT j = 3H Pe + Pj + P e + P 7 (2 24) 

dt dp e /dT 1 + dp 1 /dT 1 ' 1 ' ' 

These equations are solved simultaneously since H is a function of both T v and T7, i.e., 

p(T 7 , T v ) = p u {T u ) + p e (T 7 , B{T U )) + p 7 (T 7 ) + p B {T v ) + p 6 (T 7 ) . (2.25) 
If we define 

p = p(S = 0), e= 5 -^, x=—, (2-26) 

Po Po 
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we can estimate the effect of 9 and x on the time-temperature relation away from e + e 
annihilation: 

T 7 ~ 10 9 ^(l + 6» + X )" 1/4 t" 1/2 • (2-27) 

£ = 4.7 for three types of neutrinos. If there is no magnetic field, Eq.(2.27) reduces to the 
formula in standard BBNt 12 ! 

The modified time-temperature relationship (2.27) suggests that the contributions of 
the primordial magnetic field from both the field energy density and the electron energy 
density, accelerate the expansion rate of the universe and decrease the time scale over 
which BBN can occur. In particular, the neutrons will have less time to decay to protons 
than in the field-free case, which leads to an enhanced final ^ ratio and ultimately elevates 
the abundance of 4 He. 

Comparing the energy density from the electrons to that directly from magnetic fields, 
we find that for 7; < 10, the contribution from electrons is somewhat greater than that 
from the field, but it is still too small to be interesting with respect to the total energy 
density in the field free case. When the B field is stronger, 7$ > 10, the contribution from 
the magnetic field exceeds that of the electron phase space and dominates the total energy 
density. 

3. Limits on the field strength from BBN 

We considered all three effects discussed above in our numerical calculations to set a 
limit on the field strength allowed by BBN considerations. As expected, our calculations 
reveal that the abundances of the light elements are manifestly affected by strong magnetic 
fields (S^10 15 G). Although the three effects are important for fields in excess of B c , the 
dominant process for setting an upper limit on the magnetic field during BBN is the change 
in the time-temperature relation discussed in §2.3 in agreement with the results of Refs.[4] 
and [5]. 

Our numerical calculations show that for an initial magnetic field B(l MeV) ^10 15 G, 
the impact on the neutron to proton ratio from the magnetic field energy density, which 
decreases the neutron population, is more significant than the other effects. By using 
the observed abundance of 4 He, D, and 3 He, we find a constraint on the strength of a 
primordial magnetic field which is equivalent to an increase in the number of neutrino 
families. To explicitly calculate these effects, we set the neutron lifetime r n = 887 s ± 
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2 s,t 13 l the number of neutrino species N v = 3, and compute the primordial abundances 
numerically. Similarly to the bound on N„, the constraint on the magnetic field energy 
relies on the lower limit to r\ and the upper limit to 4 He. We use the D + 3 He lower bound 
on the baryon-to-photon ratio (77 > 2.5 x 10~ 10 ) and an upper limit to the 4 He abundance 
(Yp < 0.245) t 14 l and find that 7$ < 85. This implies that the allowed magnetic field at the 
end of BBN (T 7 = 0.01 MeV) is less than about 2 x 10 9 G which corresponds to a limit on 
the energy density of magnetic fields during BBN ps < 0.28p^. 

4. Conclusion 

In previous sections, we have provided a detailed analysis of the three major effects 
of a primordial magnetic field on the final abundances of the elements formed in big- 
bang nucleosynthesis. We have found that of the three major effects - (a) increased weak 
interaction rates; (b) enhanced electron densities in phase space; and (c) an increased 
expansion rate of the universe by the energy densities of magnetic field and electrons - 
the latter effect dominates over the modifications arising from the first two effects when 
Bi < 10 15 G even when the electron magnetic moment is included. We have computed 
these effects numerically and obtained a revised upper limit on the allowed strength of a 
primordial magnetic field on scales smaller than l^. Our results show that, in the framework 
of standard big bang nucleosynthesis, the maximum strength of a primordial magnetic field 
is such that ps < 0.28/9;,. 
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